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ABSTRACT 

We present an analysis of the central cusp slopes and core parameters of early-type galaxies using a large 
database of surface brightness profiles obtained from Hubble Space Telescope observations. We examine the 
relation between the central cusp slopes, core parameters, and black hole masses in early-type galaxies, in light of 
two models that attempt to explain the formation of cores and density cusps via the dynamical influence of black 
holes. Contrary to the expectations from adiabatic-growth models, we find that the cusp slopes do not steepen with 
increasing black hole mass fraction. Moreover, a comparison of kinematic black hole mass measurements with the 
masses predicted by the adiabatic models shows that they overpredict the masses by a factor of ~3. Simulations 
involving binary black hole mergers predict that both the size of the core and the central mass deficit correlate 
with the final black hole mass. These relations are qualitatively supported by the present data. 

Subject headings: black hole physics — galaxies: elliptical and lenticular, cD — galaxies: fundamental 
parameters — galaxies: structure 



1. INTRODUCTION 

High-resolution images obtained using ground-based and 
Hubble Space Telescope (HST) observations have shown that 
density cusps are ubiquitous in the centers of early-type galax- 
ies and true constant-density cores are extremely rare (Kor- 
mendy 1985; Lauer 1985; Crane et al. 1993; Lauer et al. 
1995). The central surface brightness profiles can be param- 
eterized by a single power law in galaxies where no core is 
resolved (power-law galaxies), while others require a double 
power law such that the slope within the resolved core is much 
shallower than the outer slope (core galaxies). According to the 
widely adopted criteria of Lauer et al. (1995) and Faber et al. 
(1997), the central surface brightness cusps can be described 
by I{r) oc r" 7 , with 7 < 0.3 for core galaxies and 7 > 0.5 for 
power-law galaxies. 

Interestingly, this classification based on the cusp slope re- 
flects a distinct dichotomy in the global photometric and kine- 
matic properties among early-type galaxies. Core galaxies are 
luminous (My < -21 mag) and have high velocity dispersions, 
compared to power-law galaxies which have lower luminosities 
and lower velocity dispersions. Core galaxies have boxy or pure 
elliptical isophotes and show slow rotation, while power-law 
galaxies tend to have disky isophotes and are fast rotators (Car- 
ollo et al. 1997; Faber et al. 1997; Ravindranath et al. 2001b; 
Rest et al. 2001). Differences between the two classes can also 
be seen in their X-ray properties (Pellegrini 1999). Thus, it ap- 
pears that the origin of the inner cusps is closely linked to the 
formation and subsequent evolution of these galaxies. 

Several models with and without central, massive black holes 
(BHs) have been proposed to explain the formation of cusps and 
the observed core properties. In models without BHs the cusp 
slope reflects the amount of dissipation that occurred during the 
formation process. Steep density cusps seen in low-luminosity 
elliptical galaxies can result from dissipative, gas-rich mergers 
(Faber et al. 1997). Numerical simulations of gas-rich merg- 
ers have shown that angular momentum transfer can lead to gas 
inflows, trigger starbursts, and result in the formation of dense 
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cores (e.g., Barnes & Hernquist 1991). Moreover, even a small 
amount of gas can cause the isophotes to become disky and in- 
duce the higher rotation velocities seen in power-law galaxies 
(Barnes & Hernquist 1996; Kormendy & Bender 1996; Faber et 
al. 1997). An analogous explanation for the weak cusps seen in 
luminous ellipticals involves dissipationless, gas-poor mergers 
(Bender, Burstein, & Faber 1992). Although such models suc- 
ceed in producing the weaker cusps seen in these objects, the 
core sizes tend to decrease or remain unchanged in the merger 
remnant (Farouki, Shapiro, & Duncan 1983; Barnes 1992). In 
fact, with the increased merger events that would be needed 
to form giant ellipticals, the cores become denser and more 
compact, in clear conflict with the empirical trend that giant 
ellipticals tend to have large-sized, diffuse cores. Even in the 
case of accretion of a low-mass secondary galaxy by a mas- 
sive primary galaxy, the dense core of the secondary survives 
the merger event, thereby resulting in a remnant that has high 
central density (Holley-Bockelmann & Richstone 1999, 2000; 
Merritt & Cruz 2001). Thus, models that aim to account for 
the formation of central density cusps without invoking BHs 
invariably produce cusps that are steeper than observed and fail 
to reproduce the scaling relations between galaxy luminosity 
and core size. 

On the other hand, numerical simulations that include the 
effects of central BHs have proved to be more promising in 
reproducing the low-density cores seen in luminous ellipticals 
(Ebisuzaki, Makino, & Okumura 1991; Makino & Ebisuzaki 
1996; Nakano & Makino 1999; Milosavljevic & Merritt 2001). 
Ebisuzaki et al. (1991) showed that the observed correlation 
between core size and total bulge luminosity can result from 
mergers of equal-mass galaxies containing central BHs. More 
recent studies have shown that the formation and subsequent or- 
bital decay of a BH binary can destroy the existing steep cusps 
of the merging galaxies via the interaction of the binary with 
surrounding stars. During the initial stages of merging, the BHs 
sink to the center under the influence of dynamical friction from 
stars and form a binary. The potential energy released as the bi- 
nary orbit shrinks can eject stars in its vicinity with high veloc- 
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ities, thereby weakening the initial density cusp and expanding 
the core (Quinlan & Hernquist 1997; Nakano & Makino 1999; 
Milosavljevic & Merritt 2001). The merger models predict that 
the core size of the merger remnant would show an almost lin- 
ear dependence with BH mass (Nakano & Makino 1999). Dur- 
ing each merger event the scouring action by the BH binary 
would eject a total mass of the order of the combined mass of 
the BHs (Quinlan & Hernquist 1997; Milosavljevic & Merritt 
2001). In addition to the BH mass, the overall mass deficit de- 
pends on the number of merger events, and is larger for the 
more luminous galaxies. Thus, a steeper-than-linear relation 
is predicted between the ejected mass and the final BH mass 
(Milosavljevic & Merritt 2001). 

An alternative model interprets the central cusps in ellipti- 
cal galaxies as the result of dynamical effects induced by the 
adiabatic growth of a central BH in an initially isothermal core 
(e.g., Young et al. 1978; Lauer et al. 1992; Crane et al. 1993; 
van der Marel 1999a). The growth of BHs in galaxy centers 
can occur by accretion of gas from stellar mass loss or other 
sources. Under the adiabatic assumption this growth is slow 
compared to the dynamical timescale, such that the phase-space 
distribution function of the stars remains unchanged. Adiabatic 
growth in an isothermal core induces a weak density cusp of 
the form p oc r" 3 / 2 , or slightly steeper for the case of a non- 
isothermal core (Young 1980; Quinlan, Hernquist, & Sigurds- 
son 1995). According to these models, the slope of the cusp is 
determined by the ratio of the BH mass to the mass of the initial 
core. A steeper cusp would result from a larger fractional BH 
mass. This is in contrast to the expectations of the merger mod- 
els, where a more massive BH would cause the cusp to weaken. 
The cusp slope and the cusp scale length of the surface bright- 
ness profiles derived from adiabatic-growth models is found to 
correlate well with the BH mass (Cipollina & Bertin 1994). 
Van der Marel (1999a) has shown that the adiabatic models can 
reproduce the range of inner-cusp slopes and core radii seen in 
the //Sr images and can be used to predict BH masses. 

In this paper, we use a large database of observed central 
parameters of early-type galaxies to test the viability of the 
adiabatic BH growth and binary BH merger models. High- 
resolution surface brightness profiles are currently available for 
many early-type galaxies. The discovery of a tight correlation 
between BH mass and bulge velocity dispersion (Gebhardt et 
al. 2000; Ferrarese & Merritt 2000) allows us to estimate BH 
masses accurately and efficiently. These recent developments 
provide an opportunity to test the models. Section 2 describes 
the sample and the sources from which the data were compiled. 
We present the observed correlations in § 3 and discuss their 
implications in § 4. A preliminary version of this work was 
presented in Ravindranath, Ho, & Filippenko (2001a). 

2. THE SAMPLE 

The sample chosen for the present study comprises early- 
type galaxies that have been observed using HST and for 
which their central surface brightness profiles (within a radius 
of ~10") have been parameterized in the literature using the 
"Nuker" function (Lauer et al. 1995). This function has the 
form 



Jb 



1+1^ 



( 7 -/3)/a 



(1) 



where (3 is the asymptotic slope as r — > oo, 7 is the asymptotic 



slope as r — > 0, ri, is the break radius at which the outer slope 
(3 changes to the inner slope 7, a controls the sharpness of the 
transition between the inner and outer slopes, and 4 is the sur- 
face brightness at r&. 

The sample is derived mainly from three sources. (1) Faber 
et al. (1997) give Nuker parameters for 58 nearby early-type 
galaxies that were observed using WFPC in the V band. These 
include galaxies in the Local Group, ellipticals and SOs in the 
Virgo cluster, and the bulges of a few nearby early-type spi- 
ral galaxies. (2) Rest et al. (2001) obtained /?-band images 
of 67 early-type galaxies using WFPC2 and fitted Nuker func- 
tions to the light profiles of 58 of these. Their sample com- 
prises all early-type galaxies from the Lyon/Meudon Extra- 
galactic Database with radial velocities less than 3400 km s" 1 , 
absolute V-band magnitudes less than -18.5, and Galactic lati- 
tudes greater than 20°. (3) Ravindranath et al. (2001b) used the 
Nuker function to perform two-dimensional fits of the central 
regions of 33 early-type galaxies observed with NICMOS in 
the H band. The galaxies were chosen from the Palomar survey 
of nearby galaxies, a ground-based optical spectroscopic study 
of a nearly complete sample of 486 galaxies with B T < 12.5 
mag and 6 > 0° (Filippenko & Sargent 1985; Ho, Filippenko 
& Sargent 1995). 

From these sources, we compiled a database consisting of 
galaxies for which the classification of the central surface 
brightness profile type is unambiguous. Only galaxies with in- 
ner slope 7 < 0.3, break radius ri, > 0'.'2, and a > 0.8 were 
included as core galaxies, to ensure that the core region is well 
resolved. Objects with 7 > 0.5 constitute power-law galaxies 
in the sample. We did not include in the present analysis the 
few objects that have intermediate slopes (0.3 < 7 < 0.5). We 
also excluded galaxies lacking stellar velocity dispersion mea- 
surements, which we later use to estimate BH masses (§ 4). The 
final database consists of 88 sources (38 core galaxies and 50 
power-law galaxies), the properties of which are given in Tables 
1A and IB, respectively. 

Since the above data were taken in different bands, it is im- 
portant to verify whether the central parameters vary with wave- 
length, as would be expected if there are significant differences 
in the stellar population or dust extinction. Ravindranath et al. 
(2001b) compared the Nuker parameters derived from optical 
and near-infrared images. Even though the outer slope (3 can be 
significantly different between the two bands (a consequence 
of radial color gradients), the inner slope 7 and the break radius 
rb show good overall agreement. Another concern when using 
data from various sources is to ensure that physical parame- 
ters are derived using the same distance scale. The apparent 
quantities from the literature were converted to absolute quan- 
tities using distances from Tonry et al. (2001), when available, 
which are based on measurements of /-band surface brightness 
fluctuations. The distances for the remaining galaxies were cal- 
culated from the radial velocities given in the Third Reference 
Catalogue of Bright Galaxies (de Vaucouleurs et al. 1991) us- 
ing a Hubble constant of Hq = 75 km s" 1 Mpc" 1 . 

3. COMPARISON OF DATA WITH MODEL PREDICTIONS 

Both the adiabatic-growth model and the binary BH model 
attempt to explain the formation of cores and weak cusps in 
elliptical galaxies based on the dynamical effects of BHs. We 
compare the relation between the observed quantities (core ra- 
dius, inner-cusp slope, and BH mass) with the relations pre- 
dicted from these models. The BH masses required for the anal- 
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TABLE 1a: Properties of Core Galaxies in the Sample 



Galaxy Name 


D 
(Mncl 


(7q 

(km s — ^ ) 


' e 


' o 


1 


J 


Reference 




(21 


(31 


(41 


(51 


(61 


(71 


(81 


A 1020 


260.0 


345 




2.480 


0.17 


1.07 


1 


A 2052 


140.8 


247 




2.457 


0.20 


1.19 


1 


NGC 524 


23.9 


245 


3.590 


2.214 


0.03 


0.72 


2 


NGC 720 


27.6 


237 


3.723 


2.634 


0.06 


1.03 


1 


NGC 1052 


19.4 


215 


3.540 


1.575 


0.11 


0.92 


2 


NGC 1316 


21.4 


250 


3.923 


1.630 


0.00 


0.16 


1 


NGC 1399 


19.9 


359 


3.611 


2.482 


0.07 


1.01 


1 


NGC 1400 


26.4 


261 


3.684 


1.626 


0.00 


0.09 


1 


NGC 1600 


62.5 


333 


4.158 


2.975 


0.08 


1.05 


1 


NGC 2832 


92.6 


328 


2.139 


2.611 


0.02 


0.74 


1 


NGC 2986 


30.7 


268 


3.789 


2.265 


0.18 


1.15 


3 


NGC 3348 


37.8 


237 


3.670 


2.258 


0.09 


0.95 


3 


NGC 3379 


10.5 


221 


3.253 


1.908 


0.18 


1.16 


2 


NGC 3608 


22.9 


195 


3.592 


1.492 


0.00 


0.28 


1 


NGC 3613 


29.1 


213 


3.536 


1.680 


0.04 


0.73 


3 


NGC 3640 


27.0 


184 


3.634 


1.761 


0.15 


1.03 


3 


NGC 4168 


30.4 


186 


3.826 


2.570 


0.14 


1.07 


1 


NGC 4261 


31.6 


309 


3.772 


2.395 


0.00 


0.00 


2 


NGC 4278 


16.0 


251 


3.406 


1.878 


0.02 


0.69 


2 


NGC 4291 


26.1 


284 


3.269 


1.784 


0.02 


0.58 


2 


NGC 4365 


20.4 


269 


3.752 


2.218 


0.15 


1.13 


1 


NGC 4374 


18.3 


293 


3.685 


2.327 


0.13 


1.12 


2 


NGC 4406 


17.1 


246 


3.875 


1.919 


0.00 


0.00 


2 


NGC 4472 


16.3 


303 


3.914 


2.317 


0.04 


0.95 


2 


NGC 4486 


16.0 


339 


3.906 


2.772 


0.25 


1.25 


1 


NGC 4552 


15.3 


263 


3.347 


1.684 


0.00 


0.03 


1 


NGC 4589 


21.9 


225 


3.643 


1.349 


0.11 


0.90 


2 


NGC 4636 


14.6 


209 


3.854 


2.388 


0.13 


1.12 


2 


NGC 4649 


16.8 


343 


3.778 


2.465 


0.15 


1.14 


1 


NGC 4874 


96.3 


266 


4.456 


3.089 


0.13 


1.12 


1 


NGC 4889 


86.5 


381 


4.120 


2.848 


0.05 


0.99 


1 


NGC 5077 


37.5 


260 


3.657 


2.467 


0.23 


1.19 


3 


NGC 5557 


42.8 


253 


3.774 


2.400 


0.14 


1.02 


3 


NGC 5813 


32.2 


230 


3.880 


2.091 


0.08 


0.99 


1 


NGC 5903 


33.9 


210 


3.762 


2.417 


0.13 


1.10 


3 


NGC 5982 


40.2 


250 


3.686 


1.954 


0.06 


0.85 


2 


NGC 6166 


121.3 


302 


4.526 


3.115 


0.08 


1.06 


1 


NGC 7768 


109.2 


286 


4.200 


2.325 


0.00 


0.00 


1 



Note. — Col. (1) Galaxy name. Col. (2) Surface brightness fluctuations distance from Tonry et al. 2001, when available, 
and otherwise derived from the heliocentric radial velocity and Hq — 75 km s -1 Mpc -1 . Col. (3) Central stellar velocity 
dispersion from Prugniel & Simien 1996. Col. (4) Logarithm of the effective radius in pc. Col. (5) Logarithm of the break 
radius in pc. Col. (6) Inner cusp slope of the surface brightness profile. Col. (7) Inner cusp slope of the deprojected 
luminosity density profile. Col. (8) Reference for 7 and rj. 

References. — (1) Faber et al. 1997; (2) Ravindranath et al 

ysis presented here are derived using the empirical correlation 
between BH mass and bulge stellar velocity dispersion (Geb- 
hardt et al. 2000; Ferrarese & Merritt 2000). Gebhardt et al. 
(2000) use a = a e , the projected, luminosity-weighted velocity 
dispersion measured within the effective radius r e . 
For our sample we do not have uniform measurements of a e , we 
use instead a = <to, the central velocity dispersion, since in gen- 
eral a e w (To within a scatter of ~10% (Gebhardt et al. 2000). 
The central velocity dispersions were taken from the compila- 
tion by Prugniel & Simien (1996). 

3.1. Cusp Slopes and Black Hole Masses 



2001b; (3) Rest et al. 2001. 

The adiabatic-growth model predicts that more massive BHs 
(relative to galaxy core mass) produce steeper density cusps. 
This is well illustrated when the inner-cusp slope (7 as param- 
eterized in the Nuker function) is plotted against the dimen- 
sionless BH mass /j, = M Bli /M mre (Cipollina & Bertin 1994; 
van der Marel 1999a), where M core is the mass of the initial 
isothermal core. We do not have an independent estimate of 
M core , since the observed core properties, by assumption, have 
been altered by the BH. However, if the initial state of the 
galaxy can be approximated by an isothermal sphere (as is done 
in the adiabatic models), the initial core radius would be a fixed 
fraction of the initial half-mass radius. Since the 
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TABLE 1b : Properties of Power-law Galaxies in the Sample 



Galaxy Name 


u 


CO 

(km s - 1 ) 


log r e 


7 


7 


.Reference 


I 1 ) 


( 2 ) 




( 4 J 


V s ) 


(6J 




NGC 221 


0.8 


77 


2.175 


0.50 


1.50 


2 


NGC 596 


21.7 


164 


3.499 


0.55 


1.55 


1 


NGC 821 


24.1 


209 


3.724 


0.64 


1.64 


2 


NGC 1023 


11.4 


212 


2.611 


0.78 


1.78 


1 


NGC 1172 


21.4 


117 


3.603 


1.01 


2.01 


1 


NGC 1426 


24.0 


153 


3.482 


0.85 


1.85 


1 


NGC 2549 


12.6 


154 


2.794 


0.67 


1.67 


3 


NGC 2592 


25.7 


282 




0.92 


1.92 


3 


NGC 2634 


33.4 


171 


3.646 


0.81 


1.81 


3 


NGC 2685 


11.7 


103 


2.955 


0.73 


1.73 


2 


NGC 2699 


26.9 


131 




0.84 


1.84 


3 


NGC 2778 


22.9 


166 


3.262 


0.94 


1.94 


3 


NGC 2872 


43.2 


271 


3.638 


1.01 


2.01 


3 


NGC 2950 


14.9 


182 




0.82 


1.82 


3 


NGC 3065 


26.6 


164 




0.79 


1.79 


3 


NGC 3078 


35.1 


237 


3.608 


0.95 


1.95 


3 


NGC 3115 


9.6 


274 


3.224 


0.73 


1.73 


2 


NGC 3384 


11.5 


164 


2.860 


0.64 


1.64 


2 


NGC 3414 


25.2 


250 


2.262 


0.83 


1.83 


3 


NGC 3599 


20.3 


79 


3.470 


0.79 


1.79 


1 


NGC 3605 


20.7 


100 


3.238 


0.67 


1.67 


1 


NGC 3900 


23.9 


115 


3.759 


0.51 


1.78 


2 


NGC 4026 


13.6 


197 


2.584 


0.51 


1.71 


2 


NGC 4121 


29.3 


118 


2.829 


0.85 


1.85 


3 


NGC 4128 


31.0 


206 




0.71 


1.71 


3 


NGC 4143 


15.9 


270 


1.790 


0.59 


1.59 


2 


NGC 4150 


13.7 


123 


1.725 


0.64 


1.64 


2 


NGC 4239 


12.5 


63 


2.989 


0.65 


1.65 


1 


NGC 4387 


21.3 


112 


3.201 


0.72 


1.72 


1 


NGC 4417 


11.2 


82 


2.134 


0.71 


1.71 


2 


NGC 4434 


26.8 


122 


3.380 


0.70 


1.70 


1 


NGC 4464 


16.5 


129 


2.632 


0.88 


1.88 


1 


NGC 4467 


19.0 


81 


2.951 


0.98 


1.98 


1 


NGC 4503 


17.9 


115 




0.64 


1.64 


3 


NGC 4551 


17.3 


114 


3.170 


0.80 


1.80 


1 


NGC 4570 


23.3 


190 


2.132 


0.85 


1.85 


1 


NGC 4621 


18.2 


237 


3.612 


0.50 


1.76 


1 


NGC 4648 


19.6 


222 


3.166 


0.92 


1.92 


3 


NGC 4697 


11.7 


174 


3.630 


0.74 


1.74 


1 


NGC 4742 


15.5 


102 


2.943 


1.09 


2.09 


1 


NGC 5017 


33.9 


174 


3.122 


1.12 


2.12 


3 


NGC 5308 


27.2 


230 




0.82 


1.87 


3 


NGC 5370 


40.7 


133 




0.62 


1.62 


3 


NGC 5812 


26.9 


213 


3.492 


0.59 


1.59 


1 


NGC 5838 


18.1 


290 


3.237 


0.93 


1.93 


2 


NGC 5845 


25.9 


241 


2.705 


0.51 


1.50 


1 


NGC 6278 


37.2 


157 




0.55 


1.55 


3 


NGC 6340 


15.9 


147 


2.380 


0.59 


1.58 


2 


NGC 7332 


23.0 


136 


2.251 


0.90 


1.90 


1 


NGC 7743 


20.7 


84 


3.321 


0.50 


1.48 


2 



Note. — Col. (1) Galaxy name. Col. (2) Surface brightness fluctuations distance from Tonry et al. 2001, when available, 
and otherwise derived from the heliocentric radial velocity and Ho = 75 km s _1 Mpc -1 . Col. (3) Central stellar velocity 
dispersion from Prugniel & Simien 1996. Col. (4) Logarithm of the effective radius in pc. Col. (5) Inner cusp slope of 
the surface brightness profile. Col. (6) Inner cusp slope of the deprojected luminosity density profile. Col. (7) Reference 
for 7 and rj. 

References.— (1) Faber et al. 1997; (2) Ravindranath et al. 2001b; (3) Rest et al. 2001. 
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FIG. 1. — The radius of the initial isothermal core for the best-fitting 
adiabatic-growth model versus the effective radius. Core galaxies are shown 
as solid points and power-law galaxies as open circles. The solid line has a 
slope of unity and corresponds to r CO re = 0.04 r e . Only galaxies common to 
this study and that of van der Marel (1999a) are shown. A representative error 
bar for r e is given in the lower-right corner of the plot. 

adiabatic growth of the BH should not significantly alter the 
mass distribution at large radii, it is safe to assume that the half- 
mass radius, or a closely related measure such as the effective 
radius, remains unaffected. In Figure 1, we plot the core radius 
from the best-fitting adiabatic model of van der Marel (1999a) 
versus the effective radius for all galaxies common to both stud- 
ies. The radius of the isothermal core indeed is well correlated 
with r e and can be approximated as r core w 0.04 r e \ the mean 
scatter about this relation is 0.23 dex. Thus, analogous to p, we 
define the dimensionless mass p! = M B n/M e , where M e is the 
mass within r e calculated using M e = 3<r 2 r e /G. The effective 
radii were taken from Faber et al. (1989), Bender et al. (1992), 
and Baggett, Baggett, & Anderson (1998). 

A plot of the observed inner-cusp slope 7 versus pi (Fig. 2) 
does not appear to be consistent with the expectation from the 
adiabatic-growth model. We have overplotted the theoretical re- 
lation between 7 and p, as given by Cipollina & Bertin (1994), 
but shifted arbitrarily in the horizontal direction to allow for a 
constant scaling between p and p! (i.e., between M core and M e ). 
The predicted line matches the data poorly. The largest discrep- 
ancy lies with the power-law galaxies, which, because of their 
steep cusps, are required to have values of p! up to 1-2 orders 
of magnitude larger than observed. Except for a larger spread, 
most power-law galaxies, in fact, are characterized by values 
of p' within the fairly restricted range found in core galaxies 
(log// « -2.9 ±0.5). In other words, the BH mass constitutes 
a roughly constant fraction of the central stellar mass within 
r e , albeit with substantial scatter. This is just the familiar lin- 
ear correlation between BH mass and bulge mass (Kormendy & 
Richstone 1995; Magorrian et al. 1998; van der Marel 1999b; 
Ho 1999; Kormendy & Gebhardt 2001). The distribution of 
points for the core galaxies might be marginally consistent with 
the model, but the scatter is large. 

Since M e is distance dependent, this test may be affected, at 
least in part, by uncertainties in the distances used. The param- 
eter 7 itself is also vulnerable to distance effects because the 
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FIG. 2. — The inner-cusp slope 7 versus the dimensionless BH mass de- 
fined as /j.' = Mftu/M e , where M e is the total mass within the effective radius. 
Core galaxies are shown as solid points and power-law galaxies as open cir- 
cles. Crosses denote galaxies for which the effective radii were derived from 
bulge-disk decomposition (from Baggett et al. 1998). Representative error 
bars are given in the lower-right corner of the plot. The solid line gives the 
theoretical prediction from Cipollina & Bertin (1994) for the relation between 
7 and fi = M-Qu/Mcom, shifted arbitrarily along the horizontal axis to allow 
for a constant scaling between ^ and //. 

cores of distant galaxies may not be sufficiently well resolved 
to yield the asymptotic inner slope of the Nuker function. This 
latter effect, however, cannot be strong for the core galaxies 
in our sample because we have explicitly chosen objects with 
well-resolved break radii (see § 2). We do not believe that dis- 
tance errors could have erased an intrinsic trend in Figure 2. 
The nominal error bars shown on the lower-right corner of the 
plot are significantly smaller than the observed spread in the 
data points. 

Figure 3a compares the BH masses inferred from adiabatic 
models by van der Marel (1999a) with independent, direct mea- 
surements obtained from stellar and gaseous kinematics (as 
summarized in Kormendy & Gebhardt 2001); a similar compar- 
ison is made in Figure 3b for a larger sample with independent 
masses derived using the M Bli - a relation of Gebhardt et al. 
(2000). The adiabatic models tend to yield masses that are sys- 
tematically larger, on average by ~0.5 dex; the offset is slightly 
larger for power-law galaxies than core galaxies. 

3.2. Core Properties and Black Hole Binaries 

Simulations of merging galaxies with central BHs produce 
merger remnants whose core properties correlate well with the 
final BH mass (Ebisuzaki et al. 1991; Nakano & Makino 1999; 
Milosavljevic & Merritt 2001). Nakano & Makino (1999) stud- 
ied the dynamical reaction of the central regions of a galaxy, 
modeled as an isothermal core, to the infall of a massive BH. 
They found that larger BHs tend to create larger and more mas- 
sive cores. Figure 4a shows the relation between the observed 
core size — here represented by the break radius r\, — and the 
BH mass. Albeit with substantial scatter, the observed trend 
qualitatively agrees with the simulations in that more massive 
BHs produce larger cores. The correlation revealed by our sam- 
ple is highly statistically significant: the Spearman's p 
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FIG. 3. — Comparison of BH masses obtained from the adiabatic-growth models of van der Marel (1999a) with {a) BH masses derived directly using stellar and 
gas kinematics (Kormendy & Gebhardt 2001) and (b) BH masses derived using the Mbh~o" relation from Gebhardt et al. (2000). The solid points correspond to 
core galaxies and open circles to power-law galaxies. Representative error bars are given in the lower-right corner of the plot. 



correlation coefficient is 0.57, with a probability of 5 x 10 -4 
that the two variables are uncorrelated. However, the observed 
slope is substantially steeper than the unity slope suggested by 
Nakano & Makino (1999). The best-fitting linear regression 
line, obtained by adopting Mbh as the independent variable 
and including nominal errors on both variables, has a slope of 
~ 2.8 ± 0.84. It would be fruitful in the future to reexamine 
this issue using simulations that incorporate more realistic ini- 
tial conditions. 

According to the merger models, more luminous galax- 
ies should have experienced a larger number of mergers than 
fainter galaxies, and hence their centers should have been 
"scoured" more thoroughly by the orbital decay of BH bina- 
ries. This scenario provides an attractive explanation for the 
observed correlation between core size and total bulge lumi- 
nosity (Faber et al. 1997; Ravindranath et al. 2001b). Recently, 
Milosavljevic & Merritt (2001) carried out more realistic simu- 
lations that start with galaxies containing central BHs and steep 
density cusps (p oc r" 2 ), as observed in power-law galaxies. Un- 
like earlier simulations where the dynamical effects of the mas- 
sive BH were studied by introducing a naked BH into an exist- 
ing core, their model includes the effect of stars bound to the 
BHs. This is found to reduce the timescale over which the two 
BHs approach each other due to dynamical friction and form a 
binary. As the binary orbit shrinks, the gravitational slingshot 
effect ejects stars from the center and causes the density cusp to 
flatten to p ~ r" 7 , where 7' < 2. Assuming that the initial cen- 
tral density profiles are represented by steep p ~ r~ 2 cusps, the 
ejected mass or the mass deficit is defined as the mass needed 
to recover the initial profile from the observed flat profile with 
slope 7'. A key prediction from the simulations by Milosavlje- 
vic & Merritt (2001) is that the mass ejected during each merger 
event by the BH binary should be comparable to their combined 
mass. The total mass ejected depends on the final BH mass and 



the number of mergers, and hence a steeper-than-linear relation 
is expected between the mass deficit and BH mass. Milosavl- 
jevic & Merritt (2001) give the following expression for the 
ejected mass: 



2(2 -V) a 2 r b 
3-7' G ' 



(2) 



where r h is the core radius (equivalent to the break radius in 
the Nuker function) and 7' is the inner-cusp slope of the lumi- 
nosity density profile. We obtained the luminosity density pro- 
file by inverting the surface brightness profile, parameterized 
by the Nuker function, using the Abel integral equation and as- 
suming spherical symmetry (Binney & Tremaine 1987). A plot 
of M e j versus M SH (Fig. Ab) shows that the two quantities are 
indeed well correlated. The best-fitting linear regression line, 
again obtained by adopting M Bll as the independent variable 
and including nominal errors on both variables, gives a slope of 
~ 2.3 ±0.51. Milosavljevic & Merritt (2001) presented a sim- 
ilar result based on a smaller sample. Formally, the correlation 
in our sample again is statistically very significant; the Spear- 
man's p correlation coefficient is 0.77, and the probability that 
the correlation could have arisen by chance is < 10" 4 . 

We stress that the tightness of the Mbh-ct relation (Gebhardt 
et al. 2000; Ferrarese & Merritt 2000) implies that there should 
be little difference between BH masses estimated in this fashion 
compared to directly measured masses. We illustrate this point 
in Figure 4 by replotting in open symbols the six objects having 
kinematically determined BH masses. 

4. DISCUSSION AND SUMMARY 

As discussed by Faber et al. (1997) and Ravindranath et al. 
(2001b), the cores of luminous elliptical galaxies populate a 
core "fundamental plane" (defined by r b , surface brightness p b 
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at rb, and cto), analogous to the fundamental plane on global 
scales (defined by r e , the surface brightness [i e at r e , and a e ). 
More luminous galaxies with higher central velocity disper- 
sions have characteristically larger, more diffuse cores. Most 

theories for the formation of early-type galaxies involve dissi- 
pative gas inflows that tend to increase the central density and 
form steep nuclear cusps. While low-luminosity galaxies with 
power-law profiles can be explained naturally in this manner, 
the existence of diffuse cores with weak cusps in massive galax- 
ies poses a serious observational challenge, as discussed by 
Faber et al. (1997). How did they form? And once created, how 
did they manage to avoid acquiring high-density cores from the 
accretion of low-mass satellites? 

A promising framework to address the formation of cores in- 
vokes the dynamical effects of merging central massive BHs. 
This is all the more appealing because there is now increasing 
evidence that BHs are ubiquitous in bulges (Magorrian et al. 
1998; Ho 1999; Kormendy & Gebhardt 2001). If so, the hierar- 
chical assembly of massive elliptical galaxies would inevitably 
lead to BH mergers. The most recent simulations by Milosavl- 
jevic & Merritt (2001) incorporate realistic initial conditions 
for the merging galaxies. Starting with two low-mass galaxies, 
each containing a central BH and a steep nuclear cusp, the hard- 
ening of the BH binary leads to the formation of a remnant with 
a large, diffuse core. We show in this paper that galaxy cores 
obey certain trends that are expected for the merger remnants 
of binary BHs. 

If cores form as an aftermath of galaxy mergers containing 
BHs, the work of Nakano & Makino (1999) indicates that the 
sizes of the cores should scale with the final BH mass. Us- 
ing the M B h - cr relation to estimate BH mass, we find that this 



indeed seems to hold in our sample: the core size increases 
roughly as r/, oc M BH (Fig. 4a). Although the slope of the cor- 
relation is steeper than predicted by Nakano & Makino (1999), 
the discrepancy should not be taken too literally given the sim- 
plified initial conditions used in their numerical calculations. 
Combining this functional dependence with Mbh oc cr 3 8 (Geb- 
hardt et al. 2000), we recover within the errors, the empiri- 
cal scaling relation between core size and velocity dispersion, 
r/, oc cr 6 (Faber et al. 1997; Ravindranath et al. 2001b). Thus, in 
this interpretation, the scaling relation observed between r/, and 
a is an indirect manifestation of two more fundamental correla- 
tions, namely the j"/,-Mbh relation which reflects the dynamical 
effects of binary BHs during core formation, and the Mbh _ cr 
relation which bears on processes associated with the joint for- 
mation of BHs and bulges (e.g., Kauffmann & Haehnelt 2000; 
Burkert& Silk 2001). 

Using the numerical models of Quinlan & Hernquist (1997), 
Faber et al. (1997) already demonstrated that BH mergers pro- 
duce cores whose masses scale approximately with the final 
BH mass. They suggested that the relation between core and 
BH mass might depend on the slope of the nuclear cusp and 
pointed out the need for simulations with more realistic initial 
conditions. The calculations of Milosavljevic & Merritt (2001) 
add a new degree of realism to the simulations by modeling 
the merging galaxies as systems with initially steep cusps, as 
seen in power-law galaxies. Adopting their formalism to com- 
pute the amount of stellar mass ejected by the BH binary during 
core formation, and again the M B h - cr relation to estimate BH 
mass, we find that M e j ex M| H (Fig. 4b). Unlike the case of 
the relation between rt, and Mbh, the slope of the correlation 
between M e j and M B h cannot be easily traced to known central- 
parameter scaling relations. The r\,-a and M m -cr relations, 
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when combined with Equation 2, yields a slope of ^4 for the 
M e j -Mbh relation, whereas the actual observed slope is ~2. 
Although each merger produces M e j w M BH , the total ejected 
mass depends on the number of mergers; more massive galax- 
ies with more massive BHs require a larger number of merg- 
ers, and thus the slope of the M e j -M Bli is expected to be > 1. 
This is in qualitative agreement with the observations. In de- 
tail, of course, the exact slope should depend on complications 
not treated in the models. These include heterogeneous mergers 
and additional processes such as gas accretion. 

A competing view of core formation posits that the nuclear 
cusps imprint the formation of the central BH through adia- 
batic growth. According to this picture, the slope of the cusp 
profile should correlate with the relative BH mass fraction; for 
a galaxy core of a given mass, a more massive BH induces a 
steeper cusp. This basic expectation, however, is not supported 
by the data. The strongest disagreement comes from the power- 
law galaxies, whose steep cusps would suggest that these ob- 
jects have exceptionally large BH mass fractions, contrary to 
what is actually observed. The core galaxies fare somewhat 
better. The distribution of cusp slopes versus BH mass fraction 
is not inconsistent with the sharp rise anticipated by the mod- 
els (Fig. 2), although the agreement cannot be said to be good. 
Regardless of this comparison, however, a basic shortcoming of 
the adiabatic-growth model is that it leaves the main problem of 
the formation of cores unanswered; the BH is assumed to grow 
in a preexisting core. 

As shown by van der Marel (1999a), one of the most attrac- 
tive features of the adiabatic model is that it potentially allows 
us to estimate BH masses based on photometric data alone. 
This would be a remarkably efficient alternative to the time- 
consuming spectroscopic observations needed for the conven- 
tional methods based on kinematics. In view of the difficulties 
described above, it might come as a surprise that the BH masses 
derived from the adiabatic models should differ from the ac- 
tual masses by just ^0.5 dex (Fig. 3). This is not entirely un- 



expected, however, given that the adiabatic models were fitted 
specifically to the core region of the central light profiles, whose 
size, rj (see Table 1 A), is typically only a factor of a few larger 
than the "sphere of influence" of the BH, r Bli ~ GM B h/c 2 - In 
other words, the mass scale of the initial core, set by the core 
size, is comparable to the BH mass. In § 3.1, we noted that the 
offset between the model-inferred masses and the actual masses 
tends to be systematically larger for power-law galaxies than it 
is for core galaxies (Fig. 3). This is a simple consequence of the 
fact that the adiabatic models need to invoke larger BH mass 
fractions (p) to fit power-law galaxies than in core galaxies, 
whereas empirically both types of galaxies share rather similar 
BH mass fractions. 

In summary, we find that the existing data on the central re- 
gions of early-type galaxies support the notion that the observed 
properties of the cores in luminous galaxies have been substan- 
tially molded by the dynamical influence of binary massive BHs 
during past merger events. The merger hypothesis can consis- 
tently account for the core fundamental-plane relations and the 
dependence of core size and mass on central BH mass. The 
observed cusp slopes do not appear to comply with the expec- 
tations of the adiabatic-growth model for BHs. 
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